Geometric parametrization of binary elastic collisions 
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A geometric view of the possible outcomes of elastic collisions of two massive bodies is developed 
that integrates laboratory, center of mass, and relative body frames in a single diagram. From these 
diagrams all the scattering properties of binary collisions can be obtained. The particular case of 
gravitational scattering by a moving massive object corresponds to the slingshot maneuver, and its 
maximum velocity is obtained. 
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I. INTRODUCTION 

We show how to geometrically parametrize the elastic 
collision between two bodies of arbitrary masses and ve- 
locities by using diagrams in the rest-frame of one of the 
bodies. For given masses and initial velocities the possi- 
ble solutions in two dimensions can be parametrized by 
a single angle 6 for both attractive or repulsive interac- 
tions. Although an elastic collision is an highly idealized 
approximation of real interactions, there are many appli- 
cations where this approximation is appropriate. 

The statement that two bodies collide means that for a 
very short time 6t in comparison to the ratio of character- 
istic length scales and speeds, the forces due to gravita- 
tional, electromagnetic, or any other interaction between 
the two bodies dominate any external forces in causing 
the momentum change Apj of each of the bodies. This 
condition implies that the impulse received from the ex- 
ternal forces by one body during the collision is negligible 
in comparison to the impulse contribution from the in- 
teraction forces due to the other body, which usually jus- 
tifies assuming that the total linear momentum or center 
of mass momentum of the two bodies during the colli- 
sion is a constant. This assumption is an approximation 
because the total change of momentum of the system is 
equal to the total impulse received. (The total impulse 
due to internal forces must add to zero because the forces 
on the two bodies must be instantaneously equal and op- 
posite.) This approximation becomes better as St or 
the external forces become weaker relative to the internal 
forces and is exact when there are no external forces in 
which case the center of mass momentum is a constant 
of the motion. 

When the interactions are conservative, the total me- 
chanical energy is conserved during the collision, which 
leads to the assumption that the total kinetic energy is 
conserved immediately before and after the collision, that 
is, when the internal forces are (and again become) negli- 
gible in comparison to the external forces. This assump- 
tion is also an approximation, which becomes asymptot- 
ically exact when there are no external forces. When the 
internal forces are central, the total angular momentum 
is also conserved in the collision under similar assump- 



tions. 

The usual treatment of the elastic collision of two 
bodies of masses mi and TO2 with initial velocities Vq 
and Uq invokes conservation of linear momentum and ki- 
netic energy in a one- or two-dimensional settingiSi^ii*^ 
Three-dimensional collisions are seldom addressed, but 
see Ref . la. 




FIG. 1: Laboratory and CM views of an elastic collision be- 
tween different bodies of different masses (m2 > mi) with 
initial velocities iio and vo- When drawn from the same ori- 
gin. Vein and the resulting velocities ui and vi always lie on 
a straight line, as also holds for Vcm and uq and vq. Further- 
more, the possible outcomes of vi and ui lie on circumfer- 
ences centered at Vcm passing through vo and uq. The angle 
if' is the scattering angle in the CM frame and 9 and if rep- 
resent the u- and v-scattering angles relative to the incoming 
velocity of the CM, viewed in the rest frame of the it-body. 



The view from the initial (rest) frame of one of the 
bodies is usually worked out and related to the center of 
mass (CM) view of the collision. The latter is particularly 
simple because the total linear momentum Pcm is always 
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zero in this reference frame, which means that 

Vq = Uq and Vj = Uj^, 

mi TOi 



(1) 



where v- = Vi — Vcm and u- = Ui — Vcm- Equation (0) 
means that the incoming velocities appear as cohinear 
opposing vectors in the CM frame and so do the outgo- 
ing velocities. Conservation of kinetic energy is a scalar 
equation which in the CM frame can be expressed as 



Vn = Vi 



and 



Un = Ui 



(2) 



The resulting velocity directions remain undefined, so 
additional information is necessary to completely deter- 
mine the velocities, for example, the scattering angle with 
respect to a reference direction. For a vector of given 
magnitude but unknown direction the possible outcomes 
define the points on a circumference centered on the ori- 
gin. In the CM frame these vectors, v'j^ and u'^, describe 
two circumferences whose diametrically opposed points 
represent the possible outcomes for the u-body and u- 
body velocities (see Fig.0. To return to the laboratory 
frame it is necessary to add the constant vector Vcm- 

In the laboratory frame the following procedure can be 
used to geometrically determine these outcomes. First 
notice that vi = Vq and Ui = Uq is also a possible solu- 
tion, corresponding to a missed collision. Extending vq 
and Uq from the origin defines points on two concentric 
circumferences whose center is pointed to by extending 
Vcm from the origin; these three points are in a straight 
line. These two circumferences define all the possible ve- 
locities in the laboratory frame, and once the direction 
of a resulting velocity is determined, then so is the other 
by the coUinearity of the three points (two on each cir- 
cumference plus the center). The collision trapezoid in 
velocity space referred to in many textbooks is obtained 
in Fig. ^ by joining all the arrowheads. 

These properties of the binary elastic collision are well 
known and will not be discussed further. Instead, we 
will develop an alternative geometric interpretation of 
the conservation laws and relate the laboratory view, the 
CM view, and the view of the collision from the initial 
reference frame of one of the bodies. This latter refer- 
ence frame is more practical because it is asymptotically 
coincident at i — > — oo with the non-inertial body frame 
of the relative coordinates and velocities for which the 
two-body problem for central forces is usually solved; 
that is, the frame in which the relative motion of the 
two masses will appear as that of a single reduced mass 
fj, = 17111712 / {mi + 7712), at the relative position of one of 
the masses, moving under forces pointing to a fixed total 
mass AI = mi+m2 at the frame's origin, where the other 
mass is at rest; this motion can be calculated for suffi- 
ciently well behaved forces. In particular, for a gravita- 
tional collision the result must be a hyperbola (unless the 
asymptotic provisos made in the preceding paragraphs do 
not apply, for example if the two bodies do not start or 
end sufficiently far apart, and then elliptic and parabolic 
collisions should be considered). 



II. ELASTIC COLLISION: TWO-DIMENSIONAL 
CASE WITH ONE BODY AT REST 

In a collision with a mass at rest at the origin, the 
resulting trajectories lie on a plane in which the total 
angular momentum Lq relative to the CM is orthogonal 
to both Vi and Ui (strong action- reaction law). 

In the frame where the M-body is initially at rest, 
the relevant conservation equations for a collision with 
a u-body with velocity Vq are 



mivo — mivi -I- m2Ui, 



(3a) 
(3b) 



which represent linear momentum conservation and ki- 
netic energy conservation before and after an elastic col- 
lision. An equivalent way of writing Eq. (O is 



mi 
m2 



(vo - vi) = ui 



(4a) 



— (vo - vi) • (vo + vi) = ui • ui. (4b) 
m2 

We use Eq. (|^ to replace mi(vo — vi)/m2 by ui on 
the left-hand side of Eq. (|4b|) . and Ui by mi(vo— vi)/m2 
on the right-hand side and obtain 

ui • (vo + vi) = ui • ( — (vo - vi) ). (5) 

VTO2 / 

Equation jSJl can be solved to express the unknown 
scalar product vi ■ Ui in terms of Vq • Ui as 



"ii — m2 ^ 
vi • ui = ■ Vo • ui. 

TOi + 1712 



(6) 



We take the scalar product of Ui with both sides of 
Eq. H4a|l and use Eq. lO to eliminate vi • Ui and obtain 



2mi 



-Vo • Ui = Ui ■ Ui. 

mi -I- m2 

Equation {Tj) expresses the magnitude 



ii = |ui| = Vui • Ul 



(7) 



(8) 



in terms of the unknown angle 6 that Ui makes with vo. 
We denote the outgoing direction of the it-body by 



Ui 

ue = — . 

Ul 

Then Eq. (jZj) is equivalent to 

2mi 



Ul 



mi + m2 



Vo • Ug. 



(9) 



(10) 



Apart from the unknown value of 6, the resulting u-body 
velocity must be 



Ul 



2mi 



mi -I- m2 



Vo cos{9)ug. 



(11) 
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FIG. 2: Collision diagram for different masses mi > 7712, u-body (mass m2) at rest uo ~ 0. The circumference with diameter 
d — 2Vcm is the locus of all possible outcomes ui . As chords of this circumference, ui and d — ui always subtend an angle 
7r/2. A choice of 6 determines ui and n, which is always orthogonal to ui and is restricted to a circumference of diameter 
l3 — 2vo. The outcome velocity vi is geometrically determined by n + ui — vo and lies on a circumference centered on Vcm 
with radius |vo — Vcml- In this diagram and represent the outcome velocities for the u- and u-bodies as seen from the 
CM frame, in which the incoming velocities would be — Vcm and vo — Vcm respectively. 



Equation (j4a|l can now be used to deduce an expression 
for vi : 

Vi=vo ui, (12) 

mi 

or in terms of 0, 

^ _ 2m2 

vi = Vo ■ Uo cos(6')ue. (13) 

mi + m2 

Equation (|13(l is not particularly illuminating with regard 
to its geometrical relation to Ui, so a more geometrical 
approach will be adopted in Sec. Ill Al for the determina- 
tion of Vi . 

A. Geometrical view 

Equation lO can be expressed in terms of the known 
vector Vcm- We set 

d = ^""^ vo = 2Vc.n, (14) 

mi + TO2 

and collect terms on the left-hand side using the fact that 
d ■ Ui = Ui ■ Ui, or 

(d - ui) • ui = 0, (15) 

that is, d — ill is always orthogonal to Ui. Equation H15|) 
shows that any admissible solution Ui will define a chord 



from the origin to a point on the circumference with fixed 
diameter defined by d (see Fig.|2}. 

The angle 9 G [— 7r/2,7r/2] formed by Ui and vo is the 
same as the angle between Ui and d. Thus, vlq is the 
unit vector in the direction determined by 9, and for a 
given choice of 9, 

ui = dcos{9) ue- (16) 

Given the geometric constraints on Ui , we can determine 
vi by substituting mi(vo — vi)/m2 for ui on the left- 
hand side of Eq. Ij4b|l . After rearranging terms we obtain 

ui ■ (vo +vi - ui) = 0. (17) 

Equation (|17|) means that the vector n given by 

n = vo+vi-Ui, (18) 

is always orthogonal to ui , 

ui-n = 0. (19) 

We will show now that, just like ui, the vector n is 
uniquely determined as soon as 9 is given. The expression 
for vi will then also be determined from Eq. H18|l as 

vi = ui - Vo + n. (20) 

We use the orthogonality condition expressed by 
Eq. (|19|l to cancel the right-hand side in the scalar prod- 



FIG. 3: Collision diagram for equal masses mi — m2, with the u-body (mass 7712) at rest (uq — 0). The ui and vi circumferences 
now coincide, and necessarily vi j | n _L ui , that is 6 + if = ^. 



uct of Eq. (|4a() with n, and replace Eq. H20() on its left- 
hand side to find (mi/TO2)n • (vg — vi) = 0, or 

n • (2vo - n) = 0. (21) 

That is, n and D — n (where D = 2vo) are always or- 
thogonal. Similar to Eq. IjlSfl. Eq. (|21l) means that n 
defines a chord from the origin to a point on the circum- 
ference with a fixed diameter defined by D. Because n is 
orthogonal to Ui, this condition defines a unique chord 
that makes the angle 4> = (7r/2) —9 with Vq. Its direc- 
tion defines the unit vector we, orthogonal to ug, and 
therefore 

n — 2vo cos — 9^\V0 — 2vo sm{9)yvg. (22) 
We use Eqs. (|22|) and for Ui and decompose vq into 



This radius is to be expected because vi = vq, Ui = 
is one possible result for the collision, meaning that the 
closest approach of the two bodies was too far compared 
to the range of the interaction forces. Confirmation that 
in general the possible vi define such a circumference 
results from verifying the orthogonality condition for two 
particular chords 

(vi - vo) • [vi - vo + 2(vo - Vcm)] = 0, (27) 

which holds when Eqs. (|25|l . 120|l . and H4a|l are used, 
together with the identities Ui ■ d = Ui • Ui and Ui • n = 0: 

(vi-vo)- ui hn = ur(ui-d+n) = 0. 

\ mi -t- 7712 ^ TTT-l 

(28) 



Vo = 7;o(cos(6')ue + sin(6')we) 
so that Eq. H20I) becomes 

'mi — 7712 



Vi Vq 



mi + 7712 



cos(0)ue -I- sin(0)we 



(23) 



(24) 



The possible outcomes of vi also have a geometrical locus 
defined by a circumference centered at 



mi 



2 7711 + 7712 
away from the origin, with radius 

7712 



Vo 



|V0 - Ve 



7711 -t- 7712 



Vq. 



(25) 



(26) 



B. Scattering angles in the u-body and center of 
mass frames 



The scattering angles for the collision can now be de- 
duced from the parameters 7n,i, 77i2, and 9. Relative to 
the invariant direction Vcm of the CM, the scattering an- 
gle for ui is evidently 9 itself. As for the scattering angle 
ip between vi and Vcm||vo, it can be calculated from 
Eqs. ^ and ^ 



cos{(f) 



Vi ■ Vq 

ViVq 



2m2 



mi+mi 



cos2(6i) 



a/i - r^^2iiri2 cos2(0)' 



(29) 



The total scattering angle between Ui and Vi is 6* -I- if. 



FIG. 4: Collision diagram for different masses mi < m2, u-body (mass 7712) at rest, uo = 0. The vi circumference is greater 
than the ui circumference so back-scattering occurs when \ip\ > |-. 



If mi = 1712, the identity H29() reduces to cos{(p) 
s'm{6), which means that (p = cf) — {Tr/2) ~ 6 and 



Ui • vi = uivi cos{ip + 9) = Q. 



(30) 



Thus, if mi = m2, the resulting velocities Ui and vi 
are orthogonal and are both chords of the circumference 
with diameter d = vq (see Fig. O. This orthogonality is 
geometrically visible from Eq. 10, which in this equal- 
mass case reduces to Vi ■ Ui = 0, indicating that the 
solutions to the collision must remain orthogonal. 

Note that if the angle ip is specified instead of the angle 
6, the determination of Ui is geometrically unique when 
mi < m2, but there is an ambiguity when mi > m2 
because there are two different magnitudes for vi with 
the same p, hence two different angles 9. Direct inversion 
of Eq. ^ yields 



cos^(^) = l 



1 + — sin^ (p) ± cos((/3) 1 
m2 ' 



-ism^((^) 



(31) 

If mi < m2 the solution with the — sign is the correct 
one. This solution can be argued by looking at the lim- 
iting case of a missed collision (p = 0) with the m2 body 
at rest. In such a case 9 is ill-defined (because the m2 
velocity remains zero), but we can see that in neighbor- 
ing collision cases, the limit oi 9 as p ^ ±0 is ±7r/2, not 
zero (see Fig.QJ. Actually = for a head-on collision, 
and because nii < m2 the outcome for Vi corresponds to 
a back-scattering with p = tt. 

If mi > m2, real solutions exist only for p e [—pL, ^l] 
(indicating that there is no back-scattering in these 



cases), where 



PL = ± arcsin 



/m2\ 
Vmi / 



(32) 



The two solutions 9±(p) obtained in Eq. H31|l now apply. 
Geometrically this result could be obtained by referring 
to Fig. Inland noting that the limiting values Vi^, p^ for 
vi and p are obtained when vi is tangent to its locus 
circumference, that is, perpendicular to v'j^ = vi — Vcm. 
Because vi^ • v'^ = 0, 



(33) 



Because |vi^| = - rl ^ (mi - 1712)/ {mi + m2)vQ, 



we 



have cos{pl) = — {ra2/miY 



or 



sin((^L) = ± 



m2 
mi 



(34) 



The scattering angle p for the w-body can also be re- 
lated to the scattering angle p' as seen from the CM- 
frame. We use Vi = v( -|- together with Eqs. 
and H26|l and obtain 



i&u{p) = 



v[ sm{p') 



7712 sin{p') 



Vcm + cos((p') mi + m2 COs(iy9') 
For equal masses 



tan((^) = tan (^-^^ ■ 



(35) 



(36) 



Likewise, the scattering angle 9 for the u-body can also 
be related to p' . If we use the fact that u[ = Vcm and 



6 



Ui = u'l + Vein, we find 

u[ sin((^') 



tan(6') 



sin((/j') 



Vcm - U'l C0S(<^') 1 - C0s((y9') ' 

which simpHfies to 

tan(6') = cot(y ), 
which is independent of the mass ratio. 



(37) 



(38) 



III. ELASTIC COLLISION: GENERAL 
TWO-DIMENSIONAL CASE 

If both bodies are initially moving when viewed from 
a laboratory frame, the same analysis can be carried out 
(see Fig. IS)). All that is necessary is to switch temporarily 
to an equivalent inertial frame Sug moving with the initial 
M-body velocity Uq. In this frame the situation is exactly 
as before, that is, an elastic collision with a body initially 
at rest. The same equations and conclusions are valid 
except that everywhere we need to let — > v,; — Uq, 
Ui -> Ui - Uq, and V, 



cm ' * cm 



Uq. 



In this case the conservation equations are 
mivo + "T-2U0 = mivi + TO2U1; 

which can be rewritten as 

mi ^ 

(vo - vi) = ui - uo; 

7712 



(39a) 
(39b) 

(40a) 



mi 

TO2 



(vq - vi) • (vo + vi) (ui - Uo) • (ui + Uo)- (40b) 



Manipulation of Eq. (|4U|I in a manner similar to Sec. |n] 
will generate the equivalent relations in this new frame 
Sua ■ Equation |(7J| now becomes 



|Ui-Uo| = 



We set as before 



2m 



mi + m2 



- — (vo - Uo) • (ui - Uo). (41) 



d = 2(V„n - Uo) = 



2mi 



mi + m2 



ue 



Ui - Uo 



-(ui - Uo), 



-(vo-flo) (42) 
(43) 



where Ug is the direction of Ui — Uq and 9 is the angle 
between Vq — Uo and Ui — Uo. Then 



2mi 



-|vo-i3o|cos(0). (44) 



|ui - uol = d • ue 

mi + m2 
Note that Eq. (|41|l is equivalent to 

(ui - Uq) • (ui - Uo) = d • (ui - Uo), (45) 
which means that 

(ui - Uo) • (ill - Uo - d) = 0. (46) 



As with Eq. (|15() . Eq. H46|l states that Ui — Uq always 
defines a chord from the origin to a circumference of di- 
ameter d. The final expression for Ui is thus 

ui = Uo H — |vo - uo| cos(0)ue. (47) 

mi + m2 

The equivalent of the n vector in Eq. (|18|) is 

n = (vo - Uo) + vi - ui, (48) 

and its orthogonality to Ui — Uo remains 

(ui - Uo) • n = 0. (49) 

Therefore Eq. H21|l becomes 

n • (2(vo - Uo) - ri) = 0, (50) 

and n also defines a chord from the origin to a circum- 
ference of diameter D = 2(vo — Uo)- Thus geometrically 
once Ui is defined then so is n and Eq. (|48|l yields the 
final expression for vi as 

vi = n + ui - (vo - Uo)- (51) 

In this general case the only invariant direction in the 
collision is that of the center of mass. It is left as an 
exercise for the reader to derive a relation between 9, ip, 
and if' in the u-body and CM frames and the scatter- 
ing deviations from the center of mass direction of the 
resulting Ui,vi velocities in the lab frame. 

IV. CONCLUSION 

We have shown how the relative velocities in a binary 
elastic collision obey simple geometric relations even for 
arbitrary masses and initial velocities. As can be seen 
from Fig. |31 the possible final velocities Ui and vi for 
given initial conditions lie on two concentric circumfer- 
ences centered at a point in velocity space defined by 
Vcm- The points defined by Ui, vi, and Vcm from the 
origin always define a straight line. These circumferences 
have radii 

ru = |ui - Vcml = ■ |vo - uol (52a) 

mi + m2 

_, TTI2 

rv = |vi - Vcml = ■ |vo - Uol, (52b) 

mi -I- m2 

where vi — Ui was replaced by vo — Uo because, from 
Eq. I|51() and the orthogonality condition (|50() (see Fig.^, 



|vo-UoP = |vo-UoP- 



|vi -uip = n- (n-2(vo-Uo))- 

(53) 

The maximum radius for either of these circumferences 
is ri = |vo ~ Uol, which occurs when the respective i- 
body mass is much smaller than the other, in which case 
the latter circumference will have a vanishing diameter, 
meaning that the velocity of the massive body is similar 
to that of the center of mass itself. If, for instance, mi » 
m2, then Vcm « vo « vi. In the limit m2/mi 0, 
— + and 

mi 



mi + m2 



|Vo - Uol 



|vo - Uol 



(54) 



FIG. 5: Collision diagram for arbitrary masses mi < m2 with initial velocities vq, uq viewed in the Laboratory frame. The 
ii-body rest frame diagram is used to obtain the final velocities ui, vi. The relation of this diagram to that represented in 
Fig. is also shown. 



All the scattering properties for binary collisions can 
be obtained from collision diagrams such as those in 
Figs. id)-©. In particular, the conditions for w-body 
back-scattering or the existence of a maximum scattering 



angle depends only on the condition r„ 



An inter- 



esting exercise would be to derive the angular 6 range for 
which an increase in outgoing velocity is obtained for a 
binary elastic collision with mi < m2, uq — 0. 

One of us has written a Java application that renders 
these collision diagrams interactively.'' The program uses 
a Live Java library developed by Martin KrauSf^ Sim- 
ulations reveal scattering situations that are not intu- 
itively obvious but can be understood when the full two- 
body motion is explored. In particular, by changing the 
asymptotic angle 9 and mass ratio milrnx, we can obtain 
the optimum incident condition for a gravity-assisted fly- 
by (or gravitational slingshot).^ It is then apparent that 
the maximum velocity attainable by the smaller mass, 
mi < m2, in a collision is 



^1, 



m2 



mi -I- m2 



|vo - uol 



(55) 



which occurs when both masses exit in the same direc- 
tion as that of Vcm- As a consequence of Eq. (|55|l . when 
the massive body is initially at rest (or when the colli- 
sion is viewed from its rest frame), the maximum velocity 
attained is ui,max = I'd- 



In future work we will show that these collision dia- 
grams are useful for calculating the eccentricities and fo- 
cal distances for open Keplerian orbits for gravitational 
scattering or repulsive Coulomb scattering. An explana- 
tion of the slingshot maneuver and gravity-assist plane- 
tary fly-by can be easily obtained. In this way the orbits 
can be viewed in the laboratory frame and a study can be 
made of the optimal incidence angle for a planetary fly- 
by that delivers the maximum velocity boost in a chosen 
direction. 
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